EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 



CERN-PH-EP-2012-237 
LHCb-PAPER-2012-015 
August 23, 2012 



Measurement of the fraction of 
Y(1S) originating from %k(lP) 
decays in pp collisions at y/s = 7 TeV 

The LHCb collaborationQ 
Abstract 

The production of XbO-P) mesons in pp collisions at a centre-of-mass energy of 7 TeV 
is studied using 32 pb^ 1 of data collected with the LHCb detector. The Xb(l-P) 
mesons are reconstructed in the decay mode XbO-P) ~^ Y(lS)j — > // + //~7. The 
fraction of T(1S) originating from Xb{^-P) decays in the T(1S) transverse momentum 
range 6 < p T r ( 151 ) < 15GeV/c and rapidity range 2.0 < y r ( 15 ) < 4.5 is measured 
to be (20.7 ± 5.7 ± 2.1~^ i ' 7 4 )%, where the first uncertainty is statistical, the second 
is systematic and the last gives the range of the result due to the unknown T(1S) 
and Xfe(l-P) polarizations. 
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1 Introduction 



The production of heavy quarkonium states at hadron colliders is a subject of experimental 
and theoretical interest [Tj. The non-relativistic QCD (NRQCD) factorization approach 
has been developed to describe the inclusive production and decay of quarkonia [2] . The 
LHCb experiment has measured the production of inclusive J/0 — > /i + fi~ [3], ip(2S) [I] 
and T(nS) — >■ ji + n~ (n = 1, 2, 3) [5] mesons in pp collisions as a function of the quarkonium 
transverse momentum p^ and rapidity y over the range < < 15GeV/c and 2.0 < 
y < 4.5. A significant fraction of the cross-section for both J/0 and T(nS) production 
is expected to be due to feed-down from higher quarkonium states. Understanding the 
size of this effect is important for the interpretation of the quarkonia cross-section and 
polarization data. A few experimental studies of hadroproduction of P-wave quarkonia 
have been reported. In the case of the Xcj states, with spin J = 0,1,2, measurements 
from the CDF jSlEj, HERA-B [8] and LHCb PHO] experiments exist, while Xbj related 
measurements have been reported by the CDF [TTj , ATLAS [T2] and DO [33] experiments. 

This paper reports studies of the inclusive production of the P-wave x&j(1P) states, 
collectively referred to as x&(lP) throughout the paper. The Xfc(lP) mesons are re- 
constructed through the radiative decay Xb(l-P) — > ^(1507 in the T{1S) rapidity and 
transverse momentum range 2.0 < y T< ^ ls ^ < 4.5 and 6 < p T r ( ls ') < 15 GeV/c. The Xb2 an d 
Xbi states differ in mass by 20 MeV/c 2 and the \bi an d Xbo states by 33 MeV/c 2 [T4]. Since 
these differences are comparable with the experimental resolution, the total fraction of 
T(1S) originating from Xb(lP) decays is reported. The results presented here use a data 
sample collected at the LHC with the LHCb detector at a centre-of-mass energy of 7TeV 
and correspond to an integrated luminosity of 32 pb" 1 . 

2 LHCb detector 

The LHCb detector [T5] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The combined tracking system has 
a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and an impact parameter resolution of 20 \mi for tracks with high transverse momentum 
(Pt)- Charged hadrons are identified using two ring-imaging Cherenkov detectors. Pho- 
ton, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers. 

The trigger consists of a hardware stage followed by a software stage which applies a 
full event reconstruction. The trigger used for this analysis selects a pair of oppositely- 
charged muon candidates, where either one of the muons has a > 1-8 GeV/c or one of 
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the pair has a pr > 0.56 GeV/c and the other has a pr > 0.48 GeV/c. The invariant mass 
of the pair is required to be greater than 2.9GeV/c 2 . The photons are not used in the 
trigger decision. 

For the simulation, pp collisions are generated using Pythia 6.4 [16] with a specific 
LHCb configuration [17]. Decays of hadronic particles are described by EvtGen [Tgj 
in which final state radiation is generated using Photos [TH]. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 
toolkit [20] as described in Ref. [21]. The simulated signal events contain at least one 
T(1S) — > decay with both muons in the LHCb acceptance. In this sample of 

simulated events the fraction of T(1S) mesons produced in x b (lP) decays is 47% and 
both the Xb(l-P) and T(1S) mesons are produced unpolarized. 



3 Event selection 

The reconstruction of the Xb(l-P) meson proceeds via the identification of an T(1S) meson 
combined with a reconstructed photon. The T(nS) candidates are formed from a pair of 
oppositely-charged tracks that are identified as muons. Each track is required to have a 
good track fit quality. The two muons are required to originate from a common vertex 
with a distance to the primary vertex less than 1 mm. 




Figure 1: Distribution of the /^ + /U~ mass for selected T(nS) candidates (black points), together 
with the result of the fit (solid blue curve), including the background (dotted blue curve) and 
the signal (dashed magenta curve) contributions. 

The invariant mass distribution of the candidates is shown in Fig. [TJ It is 
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modelled with the sum of three Crystal Ball functions [22], describing the T(1S), T(2S) 
and T(3S) signals, and an exponential function for the combinatorial background. The 
parameters of the Crystal Ball functions that describe the radiative tail of the T(1S), 
T(2S) and T(3S) mass distributions are fixed to the values a = 2 and n = 1 [5]. The 
measured T(1S) signal yield, mass and width are Nrtis) = 39 635±252, my(is) = 9449. 2± 
0.4MeV/c 2 and <Jr(is) — 51.7 ± 0.4MeV/c 2 , where the uncertainties are statistical only. 

The T (IS) candidates with a p T r ( ls ) > 6 GeV/c and a ji + /i invariant mass in the 
range 9.36 — 9.56GeV/c 2 are combined with photons to form Xb(lP) candidates. The 
photons are required to have px 7 > 0.6GeV/c and cos^* > 0, where 6* is the angle of 
the photon direction in the centre-of-mass frame of the /i + /i~7 system with respect to the 
momentum of this system in the laboratory frame. 

The Xb(l-P) signal peak observed in the distribution of the mass difference, x = 
m(/x + /i~7) — m(/i + /i~), is shown in Fig. |2] for the range 6 < p T r ( 15 ') < 15GeV/c. It 
is modelled with an empirical function given by 

= A l —^e ~^- + A 2 (x - x ) a e- {cix+C2X +C3X \ (1) 
dx V27TO- 

where Ai, AM, a, A 2 , xq, a, c\, c 2 and C3 are free parameters. The Gaussian function 
describes the signal and the second term models the background. The number of Xft(l-P) 
signal decays obtained from the fit is 201 ± 55. The mean value of the Gaussian function 
is 447 ± 4 MeV/c 2 and its width is 19.0 ± 4.2 MeV/c 2 . The expected values of the mass dif- 
ferences for the three XbjO-P) states are AM(xm) = 452MeV/c 2 , AM(xm) = 432MeV/c 2 
and AM(xbo) = 399MeV/c 2 [H]. The peak position in the data lies between AM(xb2) 
and AM(xfeo) as expected for any mixture of Xbji^P) states. 

4 Fraction of IT (IS) originating from Xb(lP) decays 

The fraction of T(1S) originating from Xb(lP) decays is determined using the following 
assumptions. Firstly, all T(1S) originating from Xb(lP) arise from the radiative decay 
Xb(lP) — > T^IS)^. Secondly, the total efficiency for T(1S) — > n + n~ as a function of 
p T r ( ls ) is the same for directly produced T(1S) and for those from feed-down from Xb(lP). 
The total efficiency includes trigger, detection, reconstruction and selection. Thirdly, the 
photon detection, reconstruction and selection are independent of the T(1S) — > 
Hence the total efficiency for Xbi^-P) is factorized as e to t(Xfc) — e cond(Xfo) ' e tot(^), where 
e tot (T) is the total efficiency for T{1S) and e con d(xb) is the conditional efficiency for Xfe(lP) 
reconstruction and selection after the T(1S) — > candidate has been selected. 

The second assumption is tested by comparing the T(1S) efficiencies obtained using 
simulated events for direct T(1S) and for T(1S) coming from decays of Xfc(l-P) states. 
These efficiencies for each pt t( - 1s ^ interval agree within the statistical error, which is less 
than 0.5%. 

The conditional Xb(lP) reconstruction and selection efficiency is estimated from sim- 
ulation as 
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Figure 2: Distribution of the mass difference m(// + /z~7) — m(fi + fj,~) for selected XbO-P) candi- 
dates (black points), together with the result of the fit (solid blue curve), including background 
(dotted blue curve) and signal (dashed magenta curve) contributions. The solid (red) histogram 
is an alternative background estimation using simulated events containing a T(1S) that does not 
originate from a Xfc(lP) decay, normalized to the data. It is used for evaluation of the systematic 
uncertainty due to the choice of fitting model. The bottom insert shows the pull distribution of 
the fit. The pull is defined as the difference between the data and fit value divided by the data 
error. 



condlW £tot(r) ^MC^) ATMC (r) ' ^ ) 

where N^(xb) and N^(T) are ^ ne number of x&(lP) and T(1S) mesons obtained from 
the fit, and N^(xb) and N^(T) are the number of generated x&(lP) and T(1S) mesons, 
respectively. The value obtained is e con d(xfc) — (9-4 ± 0.1)% for 6 < p T T ( ls > < 15 GeV/c 
and 2.0 < y T ^ < 4.5. 

The fraction of T(1S) originating from ^(lP) decays is determined from the ratio 

^prod(Xfc) _ ^V"rec(Xfe)Atot(Xfe) _ ^"rec (Xb) Acond (Xb) 

N pTod (T) ~ iV rcc (T)/e tot (T) ~ iV rec (T) ' 1 ' 

where N pro d( X b) and A r pro d(T) are the total numbers of x&(lP) — > T(lS)j and T(1S) 
mesons produced, and A r rec (x fe ) and N rec (T) are the numbers of reconstructed Xb(lP) and 
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Table 1: Number of reconstructed x&(lP) and T(1S) signal candidates, conditional efficiency 
and fraction of T(1S) originating from \bO-P) decays for different p T r ( 15 ) bins. The uncertain- 
ties are statistical only. 



p T r(1S) (GeV/c) 


6-7 


7-8 


8-10 


10-15 


6-15 


N rec {xb) 


41 ±39 


35 ±22 


91 ±30 


82 ±29 


201 ±55 


N rcc (T) 


2730 ± 64 


2193 ±57 


2866 ± 64 


2627 ±59 


10 345 ±123 


e C ond(Xfo) in % 


6.7±0.2 


8.3 ±0.2 


10.0±0.2 


12.8±0.2 


9.4±0.1 


Fraction in % 


23 ±22 


20 ±12 


32 ±10 


25 ±9 


21 ±6 



T(1S) mesons obtained from the fits to the data, respectively. As the muons from the 
T(1S) are explicitly required to trigger the event, the efficiency of the trigger cancels in this 
ratio. The fraction of T(1S) originating from Xb{lP) decays for 6 < p T r ( ls ^ < 15GeV/c 
and 2.0 < y ri - is "> < 4.5 is found to be (20.7 ± 5.7)%, where the uncertainty is statistical 
only. 

The procedure is repeated in four bins of p^ r<yis \ giving the results shown in Tabled] 
and Fig. El No significant p T r ( 1,s ) dependence is observed. The mean of the measurements 
performed in the individual bins is consistent with the measurement obtained in the whole 
p T r ( ls ) range. 

5 Systematic uncertainties 

Studies of quarkonium decays to two muons pti5|l9|IT0] show that the total efficiency 
depends on the polarization of the vector meson. The effect of the polarization has been 
studied by repeating the estimation of the efficiencies e to t(x&) and e tot (T) for the extreme 
Xb(l-P) and T(1S) polarization scenarios and taking the difference in e con d(xfo) as the 
systematic uncertainty. The largest variation is found for the cases of 100% transverse 
and longitudinal polarization of the T(1S). We assign this relative variation of 1^6% as 
the range due to the unknown polarizations. 

The systematic effect due to the unknown XbjfiP), J = 0, 1, 2 relative contributions 
is estimated by varying these fractions in the simulation in such a way that the peak 
position of the mixture is equal to the peak position observed in the data plus or minus 
its statistical uncertainty. The maximal relative variation of the result is found to be 7%. 
This value is taken as a systematic uncertainty due to the unknown Xbj(lP) mixture. 

The systematic uncertainty due to the photon reconstruction efficiency is determined 
by comparing the relative yields of the reconstructed B + — > J/ijj(K* + — > K + ir°) and 
B + — > J/ip K + decays in data and simulated events. It is assumed that the reconstruction 
efficiencies of the two photons from the 7r° are uncorrelated. The uncertainty on the 
photon reconstruction efficiency is studied as a function of p^ 1 . The largest systematic 
uncertainty is found to be 6% for photons in the range 0.6 < px 7 < 0.7GeV/c, and is 
dominated by the uncertainties of the B + branching fractions. 
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Table 2: Relative systematic uncertainties on the fraction of T(1S) originating from XbO-P) 
decays. 



Source Uncertainty (%) 

Unknown XbjO-P) mixture 7 

Photon reconstruction efficiency 6 

Signal and background description 5 

Quadratic sum of the above 10 



The systematic uncertainty due to the choice of the background fit model is estimated 
from simulated events containing an T(1S) that does not originate from the decay of a 
Xb(lP). The distribution of the mass difference obtained with these events, using the same 
reconstruction and selection as for data, is shown in Fig. El normalized to the data below 
0.38 GeV/c 2 . It describes rather well the background contribution above 0.38 GeV/c 2 , both 
in shape and level. The difference between the number of data events and the normalized 
number of simulated background events in the range 0.38 — 0.50 GeV/c 2 gives an estimate 
of the signal yield. For 6 < p T r ( ls } < 15GeV/c the signal yield obtained using this 
method is 211 to be compared with 201 ± 55 obtained from the fit. The procedure is 
repeated in each p T r ( 15 ) bin. We also study the variation of signal yield by changing 
the normalization range to 0.0 — 0.3GeV/c 2 or 0.7 — 1.0GeV/c 2 . The maximal relative 
difference of 5% is taken as the uncertainty due to the choice of the signal and background 
description. Systematic uncertainties are summarized in Table |2J 

6 Results and conclusions 

The production of %&(1P) mesons is observed using data corresponding to an integrated 
luminosity of 32pb~ collected with the LHCb detector in pp collisions at y/s = 7 TeV. 
The fraction of T(1S) originating from XbO-P) decays in the kinematic range 6 < p T r ^ < 
15 GeV/c and 2.0 < y 1 ^ 15 ) < 4.5 is measured to be 

(20.7 ± 5.7 ±2.1+ 2 ;^)%, 

where the first uncertainty is statistical, the second is systematic and the last gives the 
range of the result due to the unknown polarization of T(1S) and ^(lP) mesons. 

This result can be compared with the CDF measurement of (27.1 ± 6.9 ± 4.4)% [TT] . 
obtained in pp collisions at y/s = 1.8 TeV in the kinematic range p T r ( ls ') > 8GeV/c and 
\ V T ^\ <0.7. 

The x&(lP) decays are observed to be a significant source of T(1S) mesons in pp 
collisions. This will need to be taken into account in the interpretation of the measured 
T(1S) production cross-section and polarization. 
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Figure 3: Fraction of Y(1S) originating from %j,(1P) decays for different p T r ( 151 ) bins, assuming 
production of unpolarized T(1S) and Xb(lP) mesons, shown with solid circles. The vertical 
error bars are statistical only. The result determined for the range 6 < pt < 15GeV/c is shown 
with the horizontal solid line, its statistical uncertainty with the dash-dotted lines, and its total 
uncertainty (statistical and systematic, including that due to the unknown polarization) with 
the shaded (light blue) band. 
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